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Abstract 



■ We consider the Gaussian ensembles of random matrices and describe the normal 

' modes of the eigenvalue spectrum, i.e., the correlated fluctuations of eigenvalues 

about their most probable values. The associated normal mode spectrum is linear, 
and for large matrices, the normal modes are found to be Chebyshev polynomials of 
the second kind. We contrast this with the behaviour of a sequence of uncorrelatcd 
levels, which has a quadratic normal mode spectrum. The difference in the rigidity 
of random matrix spectra and sequences of uncorrelated levels can be attributed 
Q ' to this difference in the normal mode spectra. We illustrate this by calculating the 

£h ' number variance in the two cases. 

> 

X ' 1 Introduction 

Random matrices have been used to describe level correlations in many different areas 
of physics from nuclear levels to acoustic resonances to QCD. See, e.g., the review Q. 
The most important feature distinguishing the eigenvalues of a random matrix from a 
sequence of uncorrelated levels is the existence of strong repulsion between neighbouring 
eigenvalues. This repulsion is conveniently pictured by the Coulomb gas model adopted 
by Dyson in 1962 @. Dyson exploited the equivalence between the Gaussian ensembles 
of random matrices and a classical problem in electrostatics in which identical uni- 
formly charged parallel wires are placed on a line in a harmonic confining external field. 
The eigenvalues of the random matrix experience a repulsion which is analogous to the 
Coulomb repulsion between the line charges. This level repulsion leads to the rigidity 
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of random matrix spectra, which is most easily illustrated by the fact that the variance 
of the number of eigenvalues in an (unfolded) interval of length L grows logarithmically 
with L. This is in contrast to the linear behaviour of the number variance for a sequence 
of uncor related levels. 

All spectral properties of random matrices can be calculated from the joint pro- 
bability distribution of the eigenvalues. Given the Coulomb analogy, it is natural to 
consider the independent "normal modes" of the joint probability density, which describe 
the correlated motion of eigenvalues about their most probable values. Evidently, these 
normal modes have a simple interpretation in the Coulomb analogy where they describe 
the independent oscillations of charges on a lattice about their equilibrium positions. 
An eigenvalue can be associated with each normal mode, and these eigenvalues form 
the normal mode spectrum. In the present case, soft (hard) modes correspond to large 
(small) amplitude fluctuations in the random matrix spectrum. These collective degrees 
of freedom prove to be useful in determining the long-range spectral fluctuation measures 
of random matrices including the number variance mentioned above. 

The purpose of this paper is to calculate the normal modes of the Gaussian random 
matrix ensembles and to use them for a determination of the number variance. For 
comparison, we will also determine the corresponding normal modes and the number 
variance for a sequence of uncorrelated levels. We emphasise that the results of this 
exercise are neither new nor exact. Our intention is rather to offer a new way to regard 
the fluctuations in random matrix spectra and to emphasise the value of thinking about 
the correlated motion of eigenvalues. 



2 Normal Modes for the Gaussian Ensembles 

The Gaussian ensembles of NxN matrices have a joint probability distribution of the 
eigenvalues x\, x 2 , ■ ■ ■ , xn which is |J 

( N \ 

P N p(xi,x 2 ,...,x N ) = C N p Y[ \ x i ~ Xjf ex V \-2 N ^2 x i ) j U) 

l<i<j<N \ i=l J 

where (3 = 1,2,4 corresponds to the Gaussian orthogonal, unitary, and symplectic en- 
sembles, respectively. This distribution leads to an average level density p{x) which is 
independent of (5. In the large- N limit: 

p(x) = —V2-x 2 . (2) 
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The maxima of Pn/3 correspond to the most probable locations of the eigenvalues. It 
is sufficient to consider the single maximum with Xi < and we denote the maximum 
value of Pn/3 by P%p- One immediately obtains the following set of equations which 
determine the equilibrium positions of the eigenvalues 

Y, — ^ - Nxi = . (3) 
In the vicinity of the maximum, we approximate the logarithm of Pn/3 by 

In P N/3 = In P% p + E SxidjSxj . (4) 

The matrix C is defined as 

evaluated at the maximum. The elements of C are 



Cij — ~, ^2 • (7) 

[Xi Xj ) 

The eigenvectors of C are the normal modes of the random matrix spectrum. Each 
of these eigenvectors describes a statistically independent mode of correlated motion of 
the eigenvalues of the random matrix. Clearly, they provide the natural basis in which 
to describe the behaviour of Pn/3 m the vicinity of its maximum. 

The solution to (g) for an N x N matrix is given by the zeros of the Hermite 
polynomial, 

H N (VNxi) = . (8) 

This result was obtained by Stieltjes (see appendix A. 6 in ref. ||). It follows from the 
observation that Hermite's differential equation, 

H^{x) - 2xH' N {x) + 2NH N (x) = , (9) 

reduces to (g) at the zeros of H^. 

The results for the eigenvalues and eigenvectors of the matrix C are readily stated. 
We prove these results in sections 3 and 4. The eigenvalues of C, which satisfy 

N 

YCijSy^ = X k 8y\ k) , (10) 
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are remarkably simple. Namely, 

A fe = -kN , (11) 

where k runs from 1 to N. The i-th component of the corresponding eigenvector is 

a polynomial of order k — 1 evaluated at the equilibrium value of X{. The first four 
normalised eigenvectors can be written as 

*vP = < 12 > 

^ = Gv^if- ^ 

,(3) _ ( N-l \ 1/2 (, 2N , 



= livw^yJ l 1 -^; (14) 

^ " ^(JV-l)(JV-2)(JV-3) > / r l 2JV-3 7 ' [Lb) 

(k) 

In the large-iV limit, the eigenvector component 5y\ is given by the value of the 
Chebyshev polynomial of the second kind, Uk-i(x), evaluated at the point x = X{/y2. 
One finds that (up to corrections of order 1/iV) 

8yf = N-W U k ^ (^) . (16) 

This result is natural when one notices that the identity 

TV 
i=\ 

can be approximated for large N by 

dxp(x) 6yW(x) 8y®(x) = 5 kl (18) 

with p(x) given by (§). Eqn. (|l^) is seen to be the orthogonality relation for the Cheby- 
shev polynomials of the second kind. 

3 Eigenvalues and Eigenvectors of C for Finite N 

To prove relations ( |i~l] ) to (|i~5|), begin by looking at the definition of C stated in (|6|) and 
(0), and let C act on a power of xi. We obtain 

3=1 j& K 3 l) 
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k-2 



£(i + i)*{£; 



,k-l-2 

'3 



kx 



k-l 



£ 



Nx 



1=0 
k-2 

1=0 



1 



i# Xj 



k(k - l)x k ~ 2 - (k + l)Nxl 



Here it is understood that the sums over I are zero for k < 1 and that 



N 

E*7 



(20) 



The second identity in (1£) can be proved by induction, and eqn. (y) can be used to obtain 
the third identity. It follows from ( |l9| ) that the i-th component of the fe-th eigenvector, 

(k) 

oy\ , is a polynomial of order k — 1 in X{ and that the corresponding eigenvalue is the 
coefficient, —kN, of the highest order term, x^ -1 . 

The sums a m can be determined using the relation 

N N N 



EE r -' 

i=i j=i 



k 



8=1 



-Na k , 



(21) 



which, together with (19), lead to the recursion relation 

k-2 



Nka k = J2(l + l)cr fe _2_ 2 o"2 

1=0 



1 



k(k - l)cr fe . 



(22) 



One observes directly that ctq = N . Since the Xi come in pairs ±x (or are 0), it follows 
that &2k-i = for all k > 1. This fact and the recursion relation above permit the 
determination of all of the sums a k . 

Knowing the values of the sums, a k , it is now possible to obtain all eigenvectors of 
C using ( |l9| ) and to verify directly that the vectors (12) to (15) are indeed eigenvectors 
of C with the eigenvalues stated in (|IT|). 



4 The Large- TV Limit 

In this section we show that, in the large- iV limit, the eigenvectors of C are related to 
the Chebyshev polynomials of the second kind as indicated in equation (|l6|). Define a 
generating function G by the equation 

oo 

G(x) = o 2k x 2k , (23) 

k=0 
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where the a k are now defined by the large- iV form of equation (22) which reads 



fc-2 

Nk* k = Y,(l + l)<rk-i-2<ri • (24) 
1=0 

Multiplication by x k ~ 1 and subsequent summation over k leads to an equation which is 
equivalent to the following differential equation for G 

N^-G = (xG)^-(xG) . (25) 
ax ax 

Since o"o = N, it follows that G(0) = N, and that 

G(x) = N 1 ~ V1 ~ 2X2 . (26) 
x z 

From the definition of G in expression (^), one observes that the a 2 k are simply 

These values of a 2 k can also be obtained from the integral 

/•-v/2 

a 2k = / dx x 2k p(x) , (28) 

J-V2 

where p(x) is the Wigner semi-circle describing the average level density (^). 

With these expressions for the sums a 2 k, one can construct the eigenvectors of C 
using a Gram-Schmidt orthogonalisation procedure starting with a non-orthogonal basis 
of vectors with elements x\. The coefficients to be determined in the Gram-Schmidt 
procedure are 

-V2 



N 

a kj =-riJ"V ) = - / dx x k ~ 1 Sy^(x)p(x) . (29) 
i=i J -^ 



The structure of the integrals implies that the 5y^(x) are simply the polynomials 
orthogonal on the interval [— v2, V%] with weight p(x). These polynomials are recognised 
as the Chebyshev polynomials of the second kind, appropriately scaled. This shows that 
the elements of the eigenvectors Sy^ are given by (|l^) in the large- N limit. 

5 The Number Variance for the Gaussian Ensembles 

The small amplitude, quadratic approximation to Pn/3 m terms of its normal modes is 
useful in calculations of long range spectral fluctuation measures. We illustrate this by 
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determining the asymptotic form of the number variance for the Gaussian ensembles. 
The number variance, S 2 (L), is defined as the variance of the number of eigenvalues in 
an interval of length L. (Here, it is assumed that the spectrum has been "unfolded" so 
that the average spacing between adjacent eigenvalues is 1.) It is well-known that the 
exact number variance T&(L) for the Gaussian ensembles has the form (IJ] 

^l{L) = ^\nL + K p + 0{L- 1 ) , (30) 

where Kp is a constant. The leading logarithmic term in this expression is indicative 
of a rigid sequence of numbers. This is to be contrasted with the linear L-dependence 
which characterises the number variance for a sequence of uncorrelated levels. 

We wish to reproduce the logarithmic term in the expression for the number variance 
using the normal modes. Consider an interval of the unfolded eigenvalue spectrum 
from —L/2 to L/2. For sufficiently large N, the level density of the original spectrum 
corresponding to this part of the unfolded spectrum has the constant value of p(0) = 
v2iV/7r. Within this interval, the equilibrium position of the k-th. eigenvalue is therefore 



(0) _ 7T/C 



(31) 



k V2N 

Fluctuations in the eigenvalue spectrum will move the k-th eigenvalue to a new position 
which can be written as 

Xk = ^ + S «^i n) = ^ + ^ E ^-1 (^) • (32) 

This means that, when fluctuations are present, the eigenvalues at the (unfolded) ener- 
gies ±L/2 will have eigenvalue numbers 

L V2N „ / . ttL 



*± =± 2- T^^ra • (33) 

71=1 

Since the number of levels in the interval is now k+ — k-, it follows that the number 
variance can be approximated by the ensemble averages 

E|(L) « { ( k+ - k^) 2 ) - (k + - k.} 2 
[-] 

= ^EAK-g) (34) 

n=l ^ ' 

A L 2 - 



l -u 2 

/3vr 2 ^ n 2n ~ l \AN 

' n=l 
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Here, we have made use of the fact that terms involving Chebyshev polynomials of even 
order cancel and that the averages over the coefficients are 

(w) = S ij • ( 35 ) 
This last result can be understood by expanding the 5xi in the eigenvectors of C, 

N 

Sxi = ^2 a k5yf ] ■ (36) 
k=l 

The joint probability distribution can now be interpreted as the distribution of the 
In the approximation (||), this distribution becomes 

P N p(a u a 2 ,...,a N ) = /' v , ^ ^ ' ' ' 

= p% P n e¥Xkal > ( 37 ) 

k 

from which it is clear that the are independent and Gaussian distributed with variance 

To calculate the sum over the Chebyshev polynomials, we introduce a new variable 
9 defined as cos# = ttL/AN for L < <iN and set K = [N/2]. The terms in the sum 
can be rewritten as integrals, and after an interchange of integration and summation we 
obtain 



\ n^" 1 \W = ^9 o d9 faW • (38) 

n=l 

The value L = corresponds to 9 = ir/2. For this value, the integral equals zero. We 
can now express the approximate number variance as 



S 2 fn „ 8 [*-° Ag , S in2K9'sm2(K + l)9' 

^ {L) ~ W^~e Jo dd • (39) 



For fixed L and large K, we see that 

7T „ 7rL 



2 4(2K) 

In this limit, our approximation to the number variance can be written as 



(40) 



„,, r , 4 ri . sin 2 x 2 /7rL\ , 7r . 



where Ci(irL/2) is the cosine integral and 7 is Eulers constant. For large values of L, 
this function can be approximated by 

S^(L)«^(lnL + ln|+ 7 ) . (42) 

This expression contains a logarithmic term identical to that in (|30|), The leading 
term in the number variance for large L is thus obtained correctly by the Gaussian 
approximation to Pnb- The Gaussian approximation is not sufficient to reproduce the 
constant term in (|30|). 

6 A Sequence of Uncorr elated Levels 

It is useful to find a similar description of the normal modes for a sequence of uncorre- 
cted levels. In this case, it is easiest to proceed by considering the correlation matrix 

Dij = { (xi - (xi)) (xj - (xj)) ) = {xiXj) - {xi)(xj) , (43) 

where (x;~) is the ensemble average of eigenvalue k. This approach is slightly different 
from that adopted above when we considered a quadratic approximation to lnPjv/3- 
To the extent that this approximation is exact, the two approaches lead to identical 
eigenvectors and eigenvalues which are negative reciprocals of one another. 

A spectrum of N uncorrelated levels with unit mean level density has a Poisson 
distribution for the level spacings. The joint probability distribution for the levels can 
thus be written as a product of Poisson distributions and step functions: 

JV-l 

P N (x 1 ,x 2 ,...,x N ) = e- Xl 9{ Xl ) H e~^-^ 6(x i+1 - Xi ) 

i=l 
N-l 

= e-*» 6{ Xl ) J] 9(xi +1 - x^ . (44) 

i=l 

This form of the joint probability distribution leads immediately to the ensemble avera- 
ges (xi) and (xiXj) from which the elements of the matrix D follow: 

= min{i, j} . (45) 

(k) 

The eigenvalues, Uf., and eigenvectors, ip; , of D can be obtained by exploiting the 
special structure of D and expressing the eigenvalue problem in the following suggestive 
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manner: 







-2#) 






- 20<*> 




= u k 





(46) 
(47) 



(48) 

where 2 < i < iV — 1. Introducing the definition 0fc = (2/c — 1)tt/(2N + 1), it is readily 
verified that the normalised eigenvectors can be written as 

^ = V2N + l S[nU(Pk) ■ (49) 



1 _._-2 f ^k 



The corresponding eigenvalues are 

Wfc = ^sin" 2 ^j . (50) 

To facilitate comparison with the random matrix results obtained above, it is useful to 
multiply the io k by [p(0) ] -2 = it 2 /2N 2 to establish identical scales. We then see that the 
hardest normal mode eigenvalue for the uncorrelated levels is of order —2N 2 /(t{ 2 lun) = 
— (8/ir) 2 N 2 which is comparable to the hardest eigenvalue of Xjy = —N 2 obtained for 
the random matrix ensembles. (The difference is not significant and is largely due to the 
limitations of the quadratic approximation for the joint probability density.) Significant 
differences are found in the nature of the soft spectrum. For small values of k, the 
uncorrelated levels reveal a quadratic spectrum, —2k 2 , which is in sharp contrast to the 
linear spectrum, —kN, of the random matrix ensembles seen in expression (|Tl]). 
The number variance for a sequence of uncorrelated levels is well known [Q| 

T?(L) = L. (51) 

Here, we wish to proceed in the spirit of the previous section and reproduce this result 
using the normal models just obtained. Considering the interval [0, L] and repeating the 
arguments which led to equation (|34|), we arrive at the expression 



1 N 1 

ra ^TTT £ ■ 2T7 /9 x sin 2 (L<^) • (52) 
2N + 1^ sm A {<j> k /2) 

This sum can be performed exactly when L is an integer or a half-integer. For these 
cases, we find that the number variance is given as 

. „ for integer L 

S 2 (L) PS { (53) 
- 2 (2JV+i) f° r half-integer L 
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for physically interesting values L < N. This result agrees with ( pl|) in the large-iV 
limit. For other values of L, it is useful to approximate the sum in ( |52|) by an integral 
as in eqn. (|4l|). 

9/ , 1 /"V 2 sin 2 (2Lx) 1 r , , 

E 2 L)«-/ dx ^ = 2vrL Si 2vrL) - 1 + cos 2vrL , (54) 

where Si(2-7rL) is the sine integral. For large values of L, this leads to the approximation 

S2(L) ^ L " 2^1 (2?rL) " i ' (55) 

where we note that the final term of —1/tt 2 is an artifact of having replaced sin 2 x by 
x in the denominator of (|5~4]). Evidently, eqns. (p^|) and ( |55| ) invite comparison with 
the results of eqns. ( fH]) and (^) obtained for the Gaussian ensembles. We see that 
the quadratic approximation to the joint probability density also provides a reliable 
description of the number variance in the large L limit for uncorrelated levels. 



7 Discussion and Conclusions 

We have considered the small amplitude normal modes describing the fluctuations of 
the eigenvalues of random matrices about their equilibrium positions. In the limit of 
large matrices, these modes are essentially plane waves. (Recall that the Chebyshev 
polynomials have the form U2 n (x) = cos (2n + l)x and U2 n +i{x) = sin (2n + 2)x in the 
limit of large n and fixed x.) The mean square amplitude for each mode is inversely 
proportional to its associated eigenvalue, see eqn. (|37|). Since these eigenvalues grow 
monotonically with wave number for both the Gaussian ensembles and for uncorrelated 
levels, longer wave length fluctuations have a larger amplitude. Thus, the most probable 
fluctuation in a random matrix spectrum corresponds to a common shift of all eigenvalues 
with no change in their relative separation, see (|l2|). The next most probable fluctuation 
is a simple "breathing mode" of the spectrum, see (|l3|). 

The properties of the normal mode spectrum provide us with some insight regarding 
the qualitative behaviour of long-range spectral measures. For uncorrelated levels, the 
long wave length spectrum is particularly soft with a quadratic spectrum. As a result, 
S 2 (L) is completely dominated by soft modes when L is large. This is seen most easily 
from eqn. (|54|). The linear asymptotic behaviour of the number variance for uncorrelated 
levels is thus seen to be a direct consequence of the quadratic dispersion relation obeyed 
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by the soft modes. The situation is qualitatively different in the case of the Gaussian 
ensembles where all modes obey an exact linear dispersion relation. From eqn. (|4l| ) we 
see that this linearity ensures that S 2 (L) must grow logarithmically for large L and 
that all normal modes contribute democratically to this asymptotic behaviour. Similar 
linear dispersion relations characterise perfectly elastic solids, and it seems useful to 
regard the spectral rigidity of random matrices as a consequence of the physical rigidity 
of a classical one-dimensional array of line charges. Dyson was led to the same conclusion 
|H . Although aware of the somewhat arbitrary nature of distinctions between phases in 
one dimension, he felt it appropriate to call the Coulomb gas (and hence the spectrum 
of a random matrix) a "crystal" . The considerations presented here provide additional 
support for this designation. In the same spirit, the quadratic dispersion relation for 
the normal modes of uncorrelated levels leads us to regard them as a gas. 

We have emphasised that fluctuations in random matrix spectra are most naturally 
described as the highly correlated motion of individual eigenvalues, i.e., the normal 
modes. Given the product form of eqn. (|37|), it is also clear that these normal modes are 
statistically independent. If the matrix D of eqn. (|43| ) is calculated in numerical simula- 
tions, the statistical errors associated with its eigenvalues will become uncorrelated as 
the sample size increases. Soft modes (with large amplitudes) and their eigenvalues can 
thus be determined accurately with relative ease. This fact can be useful in numerical 
simulations of ensembles which do not readily permit analytic analysis. 
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